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SIG"ir ICANCE  AND  EXPLANATION 

We  study  the  homogenization  of  a  periodically  distributed  mixture  of 
two  (or  more)  pure  media  which  can  undergo  one  or  more  changes  of  phase. 
The  problem  is  to  find  the  properties  of  the  idealized  homogeneous  medium 
(if  it  exists)  whose  behaviour  the  original  mixture  approaches  when  the 
periodic  mesh  gets  finer  and  finer. 

It  is  proved  under  quite  general  boundary  conditions,  and  under  the 
restriction  that  for  each  original  medium  the  heat  conductivity  is  inde¬ 
pendent  of  the  temperature,  that  such  a  limit  medium  exist  and  that  it  has 
the  following  properties:  its  heat  diffusion  is  homogeneous  but  non - 
isotropic  in  general,  its  specific  heat  is  an  average  of  the  specific 
heats  of  the  original  media,  and  it  does  undergo  change  of  phase  at  all 
the  temperatures  of  change  of  phase  of  the  original  media.  This  last  fact 
is  in  strong  agreement  with  everyday  experience,  as  witnessed  for  example 
by  the  freezing  of  ice  cream. 


HOW  TO  HOMOGENIZE  A  .NONLINEAR  DIFFUSION  EQUATION 
(e.g.  STEFAN'S  PROBLEM) 

Alain  Damlamian 

1.  Introduction 

For  many  physical  studies  of  composite  materials,  one  is  more  interested  in  a  global 
or  "macroscopic"  behavior  of  a  composite  medium  rather  than  in  a  detailed  "microscopic" 
one.  To  put  it  is  a  different  way,  and  for  reasons  that  can  also  come  from  numerical 
analysis  when  discretization  is  considered,  one  is  interested  in  finding  the  relevant 
properties  (i.e.  the  constitutive  laws  or  physical  parameters)  for  an  idealized  homo¬ 
geneous  medium  which  would  have  the  limit  behavior  of  the  composite  material  when  the 
size  of  the  periodic  mesh  goes  to  zero. 

Finding  the  relevant  parameters  of  this  idealized  limit  (when  it  exists)  is  the 
origin  and  one  of  the  main  contributions  of  homogenization  theory.  Since  this  is  not 
intended  to  be  an  introduction,  let  alone  a  survey, of  homogenization  theory,  the  reader 
is  referred  to  Bensoussan-Lions-Papanicolaou  [1]  for  a  complete  set  of  references. 

The  model  problem  we  are  looking  at  here  is  the  problem  of  homogenization  of  a 
non  linear  heat  equation  for  a  composite  material  consisting  of  a  periodic  mixture  of 
media  which  can  undergo  changes  of  phase,  (Stefan's  problem).  One  can  also  apply  the 
present  results  to  electromagnetic  composite  materials  (see  Bossavit-Damlamian  [1]. 

The  plan  is  as  follows 

2.  The  model  e-problem;  weak  formulation. 

3.  Estimates  for  the  e-problem. 

4.  A  short  review  of  elliptic  homogenization. 

5.  The  limit  problem  and  its  constitutive  laws. 

Conclusions 
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to  Professor  C.  Tartar  and  F.  Murat  for  their  helpful  suggestions. 
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2.  The  Modal  e-problem. 

In  short,  the  problem  we  are  looking  at  is  the  following: 

N 

Let  n  be  a  given  bounded  domain  (in  ]R  usually  N  =  3)  with  smooth  boundary  .  . 

We  restrict  ourselves  to  two  media  and  M  •  Their  distribution  in  is  given 

according  to  a  periodic  structure  of  mesh  size  e  proportional  to  a  basic  period  Y  of 
size  1.  The  basic  period  Y  is  partitioned  into  two  smooth  subsets  Y^  and  Y^  corre¬ 
sponding  to  each  medium  and  ,  so  that  correspondingly  Q  is  partitioned  into 

..^  and  ^  .  The  boundary  between  Y^  and  Y^  is  denoted  by  1  ,  and  its  image 

in  ..  is  2^  which  is  supposed  to  be  rigid  and  perfectly  heat  conducting.  As  for 
3..  =  T  it  is  split  into  F^  ^  and  F^  ,  ^  ,  corresponding  to  each  medium. 

It  is  assumed  as  customary  for  a  Stefan  problem  that  the  variations  in  volume  are 
negligeable. 

Time  will  be  restricted  to  an  interval  [0,T]  and  it  will  be  shown  that  the  result 
is  independent  of  T  . 

In  2  =  ]  0,T (  *  .I  ,  each  change  of  phase  for  each  medium  will,  in  the  strong  formula¬ 
tion,  generate  a  free  boundary  separating  the  phases.  With  our  convention  that  there  is 

at  most  one  change  of  phase  for  each  medium,  two  boundaries  S.  (i  =  1,2)  (not 

1  /  £ 

necessarily  connected)  are  generated. 

Before  writing  the  strong  formulation,  we  introduce  some  notations:  for  each  medium 
i  (  =  1,2)  a  (v) ,  th  ,  b^,  ,  denote  the  specific  heat  (a  function  of  the  temperature  v  ), 

the  temperature  of  change  of  phase,  the  latent  heat,  the  heat  conductivity  (which  is 
assumed  to  be  strictly  positive  independent  of  the  temperature  v  ,  a  limitation  for  our 
final  result  but  which  has  not  been  lifted  so  far) . 


We  can  now  write  the  strong  formulation  of  the  pron l*-m.  Sn»-  look;;  1  or  v  (t,x) 


(temperature)  and  the  two  surfaces  S  satisfying 

i ,  e 


(2.1) 


in  (0,T)  x  nije  \  Si(C 


a.  (v  )  -r— 
i  3t 


div  (k .  7  v'  ) 
1 

where  f  is  an  internal  heating 
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on  (0,T) 


.  continuity  of  v 


.  continuity  of  heat  flux 


(k  ;  v  •  n  =  7  v  •  n  ) 


b.  cos(n.t)  -  >  [h.  ■” —  cos(n,x.)}  *  0 

i  V  i  3x.  ^  S. 

3  3  J 

where  n  is  the  unit  normal  to  S.  in  space  time,  [  ]  indicates  the 

si,, 

jump  across  S.  along  n  .  This  is  the  classical  Stefan  condition  on  the 

e 


free  boundary. 


Initial  condition:  v“(0)  given  in  S)  together  with  the  initial 

(2.4)  -j  boundaries  ^  (0 )  ;  -hey  are  assumed  to  be  compatiLle  (vL(0)  = 

on  s.  (0)  . 

I 

As  for  the  lateral  boundary  conditions,  it  is  known  in  the  linear  homogenization  theory 
that  provided  they  are  of  variational  form,  they  do  not  interfere  with  the  limiting 
process.  To  be  complete  we  shall  take  them  linear  inhomogenous  of  mixed  type.  We 
assume  a  smooth  partition  of  T  in  I’+  and  T  (T  with  non  empty  interior)  and 


require  the  following: 


(2.5)  v  (t,x)  =  g  (t,x)  on  T 

-  c 

(2.6)  k.  +  P  vC  =  g+(t,x)  on  T  n  T. 

i  on  l ,  e 

here  P  is  a  non-negative  smooth  function  measuring  the  permeability  of  the  boundary 
r+  to  heat  flow,  g  and  g+  are  given  smooth  functions.  It  is  also  assumed  that 
the  boundary  data  g  ,g+  agrees  with  the  initial  data  v^  at  t  =  0.  It  turns  out 
that  conditions  (2.2)  and  (2.3)  are  Rankine-Hugoniot  type  conditions  for  the  energy 
balance  equation  taken  in  the  distribution  sense  on  Q  .  In  order  to  write  this  equation 
(which  gives  a  weak  formulation)  we  need  some  notations: 

For  each  i  ,  let  -y\  denote  the  maximal  monotone  graph  defined  (up  to  a  constant) 


V.  (v)  =  /V  a.  (s)ds  +  b.  H(v  -  6.  ) 
l  J  i  i  l 
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(H  is  the  Heaveside  function) .  This  represents  the  enthalpy  as  a  function  of  the 
temperature. 

We  also  put: 


(2.8) 

y(y,v)  =  -^(v)  for  y  in  y 

(2.9) 

k(y)  =  ki  for  y  in  Y 

Then  (2.1), 

(2.2),  (2.3)  reduce  to 

(2.10) 

-  div  (k (*)  V  vL)  =  f  in  £>•  (Q) 

JL  £ 

(2.11) 

o'"  ( t ,  x)  f  y(  —  ,  v £  (t , x) ) 

The  initial  conditions  (2.4)  ,  can  be  expressed  in  terms  of  u  alone  as  an 
initial  condition 


(2.12)  uC(0)  =  uQ  , 

which  we  can  assume  independent  of  e  . 

For  the  lateral  boundary  conditions,  we  introduce  an  auxiliary  problem,  where  the 
time  t  is  a  mere  parameter: 

Let  g£(t,x)  be  the  solution  of 

-  div(k  (  *  )  Vgc(t))  =  0  in  fi 
gE(t)  =  g  (t)  on  T 


(2.13) 


k<  +  P  g£(t)  =  g+(t,  on  r+  . 

c  on 

Clearly  gE  is  bounded  in  uniformly  in  e  .  Now  (2.1)- (2.6)  has  the 

following  weak  formulation 

(2.14)  /  -  if'  uE  +  /T  a£(v£  -  g £,*)  =  J  i<t  +  /  tf(0)u 

Q  0  Q  a 

for  all  in  C1 (Q) ,  *  (T)  =0  and  *  =  0  on  (0,T)  x  r“ 

(2.15)  v£  -  g''  =  0  on  (0,T)  x  r”  . 

In  (2.14)  a£  is  the  bilinear  Dirichlet  form  given  by 

(2.16)  a£(w,^)  =  J  k(  —  )(7wW)  +  /  p  w  •  . 

«  E  .  +  r+ 

The  smoothness  assumptions  made  for  T  ,r  allow  for  (2.14)  to  take  £  in  a  larger 

12  1 
class  namely  y?  c  W  '  <0,TjV)  where  v  is  the  variational  space:  c  H  (fi),c/j,-  =  0). 

(see  Damlamian  [1]  for  a  detailed  study). 
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(2.17)  Definition:  (u£,v£)  is  a  weak  solution  for  problem  (2 .1)  —  (2.6)  if  and  only 

if  they  satisfy  (2.11),  (2.14),  (2.15). 

(2.18)  Theorem  (see  Damlamian  (1)  and  [2]). 

Under  the  hypothesis  that  the  a^'s  are  bounded  above  and  below  away  from  zero, 
there  exists  a  unique  solution  (u£,v£)  for  problems  (2.11),  (2.14),  (2.15)  which 
satisfies : 

u£  €  W1'2(0,T;V*)  n  l“(0,T;L2((2)) 
v£-g£  £  w1'2  (0,T,L2  (fl) )  n  l“(0,T;V). 

Instead  of  giving  a  detailed  proof  we  will  only  give  the  idea  of  how  to  obtain 
uniform  estimates  in  the  following  paragraph. 
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3.  The  uniform  estimates. 

(3.1)  Proposition.  The  solutions  (u^.v1)  e  >  0  satisfy  the  following: 

12  oo  2 

uE  is  bounded  in  W  '  (OT  ;V*)  n  L  (OT  ;L  (.1))  uniformly  in  >  0  , 

v~  -  g£  is  bounded  in  W* ' 2  (OTj L2  (:1> )  n  L  (OT,V)  uniformly  in  r  •  0  . 

Proof :  To  obtain  these  estimates,  it  is  enough  to  show  them  in  the  case  of  smooth  >. 
Then  ue  and  vL  are  smooth  enough  to  replace  (2.14)  by 

(3.2)  /  -  y'»ue  +  /T  a£(v£-  ge,*)  =  /  U  +  /  *(0)u  -  /  s?  ( T ) u  ( T ) 

Q  0  Q  o  Q 


1.2 

for  all  \fi  in  W  (0T;V). 

Then  taking  $  -  [a  ]  (A**  being  the  operator  associated  to  a£  on  V)  one  gets: 


(3.3)  tlluc(t)ir  *  +  /  /  Vcuc  <_  ill  U  It  +  / 

'on  l'v  (o,t)xf. 


C  C 

U  g  + 


/*  lluell  .llfll  . 
q  e,V*  e,V* 


where  II  ^  is  the  dual  norm  of  (a  ‘  on  V  r  the  latter  being  uniformly 

equivalent  to  the  standard  norm  on  V  . 

Also  taking  i  =  [aC]  1  one  gets 

at 


(3.4)  c  |vc(t)|2 

L  (SI)  0 


ft  du  ||  +  — 

]  E  v  C2  (a  constant  which  depends  upon  f ,VQ,g  ,g  . 


From  (3.3)  and  (3.4)  one  gets  (because  v  u  can  be  assumed  non  negative),  that  v 


is  bounded  in  L  (0,T;L  (Q) ) 

e  1,2. 

u  in  w  (0,T;V*) 


Then  one  takes  >e 
rt  idv£ 


d  (vE-gc) 
dt 


to  get 


(3.5)  c  /  |—  |  dt  +  aL  (vC-gE,vE-gC)  <  c  (a  constant  which  depends  upon 

J  0  L^(fi)  4 


f»g+,g  ,vo ' • • •  ) • 

e  12  2 

From  (3.5)  one  infers  that  v  stays  bounded  in  W  '  (0,T;L  ( Q) )  and 
ve-g£  stays  bounded  in  L°°(0,T;v). 

A  detailed  proof  of  the  above  can  be  found  in  A.  Damlamian  [1]  and  [2] ,  and  in  much 


simpler  cases  in  Brezis  (4)  and  Lions  [1] 


It  is  worth  noticing  that  given  the  estimates  of  (3.1)  (even  not  uniform  in  ■  ) 


equation  (2.14)  can  be  replaced  by  (3.2)  or  even  by 

(3.6)  /fc  at(v’-g”,v)  =  /  f*i  +  /  v  (u^  -  u'(t)) 

0  Q 

for  all  v  in  V  (independent  of  t  ) . 

This  remark  (cf.  Damlamian  [1))  shows  that  (u  ,v")  is  the  solution  of  a  simpler  varia¬ 
tional  inequality  (of  the  type  studied  by  G.  Duvaut  11)). 

Another  way  of  looking  at  (3.6)  is  the  following,  since  the  operator  does  not  depend 
upon  time 

f  a'  (Vl(t),v)  =  /  (F  (t )  +u-u  (t))v> 

(3.7)  I 

|  V*  (t)  =0  on  r 

for  all  t  •  (0,T),  all  »  in  V  ,  where 

V>_(t,x)  =  /t  (v£(s,x)  -  gt(s,x))ds 
0 

F (t ,x)  =  /t  f (s ,x) ds. 

0 
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4.  A  short  review  of  elliptic  homogenization 

The  purpose  of  this  paragraph  is  to  show  how  elliptic  homogenization  works  and  hov 
it  can  be  applied  in  the  present  problem.  See  Bensoussan-Lions-Papanicolaou  [1]  (Alsc 
L.  Tartar  [1]).  With  the  same  notations  as  above,  we  consider  the  operator 
A  -  -div(k  (  —  )V)  on  U  .  Let  w£  be  the  variational  solution  of 


G 

w  »  g  on  i 

w  x  .  3w£  g  +  + 

k(  —  )  -r—  +  p  w  *  g  on  r 

£  dn 


(4.2)  aE(w L  ,'f)  =  I  f'V  +  /  g+  <e 

n  r+ 

for  all  in  V  ,  wE  =  g  on  r 

e  0  2 

We  assume  that  f  converges  to  f  weakly  in  L  (ft) . 

(4.3)  Proposition:  As  e  goes  to  zero,  we  converges  weakly  in  H1  ( .1 )  to  the 

solution  w°  of  the  following  problem 

(4.4)  a°(w°,v>)  =  j  f°vJ  +  /  g+lf 

a  r+ 

for  all  ;  in  V  ,  w°  =  g  on  T  ,  where  a°  is  the  bilinear  form  given  by 
0  ,  .  r  r  3w  r 

a  (w,4>)  =  /  1  57*'/^ 

Q  j,t  dxj  dxl  r 

with  constant  coefficients  q.  given  by 

3  f*- 


dj,Z  mes(y) 


/  My)  v  (X*  ye)7<xj  -  y  ) 


where  xJ  is  the  solution  (defined  uniquely  up  to  a  constant)  of 

(4.6)  -div(k(y)  V  x^)  =  -div (k (y ) e , )  ,  x^  periodic  in  Y  . 

N 

e..  is  the  j’th  unit  vector  in  IR  ,  y^  being  the  coordinate  on  e.  . 

Proof :  It  is  clear  that  w£  being  bounded  in  nhft)  (by  coerciveness  of  aL  with  the 

Dirichlet  boundary  condition) .  So  we  can  assume  (via  uniqueness  of  the  solution  for  the 

limit  problem- to-be)  that  we  -  w°.  Then  by  a  result  of  Tartar  [11,  one  can  see  that 

k  (  —  )  converges  weakly  in  L?  (ft)  (hence  in  (f.) )  to  J  q.  .  (q.  .  given 

e  3x  3  loc  .  j!  5x,  je 

*■  3  3 

by  (4.5),  (4.6)),  so  that  (4.2)  goes  to  the  limit  to  (4.4),  which  is  the  weak  formulation 
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(4.7) 


,0  0  0 

A  w  =  f 


3w 


3y 


+  p  w 


37 


is  the  conormal  derivative  for  A  ,  which  one  should  notice  is  not 


diagonal,  but  still  symmetric,  and  with  constant  coefficients 

2 

0  r  3 

A  ~  ~  -K  “V  3x  .  3x  * 

],2  J  3  i 


Here  we  have  also  used  the  compactness  of  the  trace  operator  from  H  (fi)  into  L 
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5.  The  limit  problem  and  its  constitutive  laws. 


Making  use  of  the  results  of  paragraph  4,  one  sees  that  g'(t)  converges  weakly  in 


H  (s'.)  to  the  solution  g  (t)  of 


0  0,  . 

A  g  (t)  =  0 


(5.1) 


g  (t)  =  g  (t)  on  " 

3q°  (t)  0  ,  + 

JT. -  +  p  g  (t)  =  g 

31  0 

A 


Also,  by  the  estimates  of  (3.1)  one  can  extract  a  sequence  of  values  of 


(*) 


going 


to  zero  such  that 

uC  -*  u°  in  W1,2(0T,V*)  n  L°°{0T;  L2  (K) ) 
vE-gE  -  v°-g°  in  W1 ’2 (0,T,L2 (Q)  n  L°°(0T,V)  . 

Hence  u  converges  strongly  in  C([0,T];V»)  and  for  all  t  •  10, T],  u’ (t)  converges 

weakly  to  u°(t)  in  L2  (SI)  . 

Consequently  we  can  apply  the  result  of  (4.3)  to  (3.7)  so  that  VE(t)  which 
obviously  converges  to  V°(t)  =  /t  (v°(s)  -  g°(s))ds  satisfies: 


(5.2) 


°<V°<t),,)  =  /  (F(t)  +  u  -  u°(t)K 

a  0 


V  (t) 


on  r 


Using  the  equivalence  with  the  weak  formulation  of  type  (2.14)  we  get 


(5.3) 


0  0 

/  -  g  =  0  on  (0,T)  *  F 


f  0  ft  0  0  0  t  t 

j  -  <f*  u  +  J  a  (v  fV>  +  J  u  (0) 


for  all  f  in  W1,2(0T;V)  f  (T)  =  V  .  .  _  =  0  . 

1 (0T)x 1 

We  now  turn  to  (2.11),  that  is 

uE(t,x)  e  y(  |  ,  vE(t,x))  . 

Assuming  e  is  chosen  so  that  gE  converges  to  g°  weakly  in  H^fO),  then  vE 
converges  to  v°  weakly  in  W1'2  (0,T;L2  (ft)  n  L°° (OTjH1  (J2) )  so  that  the  convergence  is 


<*) 

The  limit  problem  having  a  unique  solution,  it  will  be  clear  by  the  end  of  the  proof 
that  the  whole  sequence  converges. 
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6.  Conclusion 


We  conclude  that  the  limit  equations  correspond  to  the  weak  formulation  of  the 

following  strong  problem: 

du  0 
—  +  A  v  =  f 


(6.1) 


u(t,x)  e  y(v(t,x)) 

u (0 ,x)  =  u°(a) 

v(t,x)  =  g  (t,x)  x  c  T 

3v  (t,x)  +  p  v(t,x)  =  g+(t,x)  x  f.  T+ 


3v 


or  to  look  at  it  from  the  Stefan  problem  point  of  view,  a  non  isotropic  Stefan  problem. 


One  can  notice  that  we  recover  the  heat  diffusion  operator  of  the  linear  case, 
that  is  a  homogeneous  but  anisotropic  heat  diffusion. 

One  also  gets  an  "averaging”  phenomenon  for  the  graphs  y. 's  over  Y  ,  which  is 
the  only  averaging  consistent  with  the  fact  that  both  y^'s  are  defined  up  to  an  additive 
constant  and  so  is  y  .  Both  temperatures  of  change  of  phases  appear  for  discontinuities 
of  y  ,  which  is  in  agreement  with  daily  experience  (any  other  averaging  of  Y^  and  y^ 
would  have  yielded  no  discontinuity  in  the  average,  hence  no  change  of  phase).  It  is 
easy  to  see  that  the  specific  heat  and  latent  heat  of  the  limit  medium  are  averages  over 
Y  of  the  corresponding  terms . 

Finally,  on  the  theoretical  side  of  things,  it  is  of  interest  to  realize  that  the 
isotropic  diffusion  laws  are  not  stable  under  homogenization  of  Stefan  problems,  but 
anisotropic  ones  are  stable. 

It  remains  to  prove  that  the  above  can  be  extended  to  the  case  of  temperature  — 
dependent  heat  conductivity  for  each  medium,  which  is  already  more  complicated  but  solved 
in  the  non  Stefan  case  (see  Tartar  (!]). 
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